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Abstract 



We conjecture that in one-dimensional spatially extended systems the 
propagation velocity of correlations coincides with a zero of the convective 
Lyapunov spectrum. This conjecture is successfully tested in three different 
contexts: (i) a Hamiltonian system (a Fermi-Pasta-Ulam chain of oscillators) ; 
(ii) a general model for spatio-temporal chaos (the complex Ginzburg-Landau 
equation); (iii) experimental data taken from & CO2 laser with delayed feed- 
back. In the last case, the convective Lyapunov exponent is determined di- 
rectly from the experimental data. 
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Since the recognition that deterministic chaos is an ubiquitous feature of nonhnear sys- 
tems, much has been understood about their dynamical properties. An important example 
is the discovery of the relationships between Lyapunov exponents, measuring the divergence 
rate of nearby trajectories, and either geometrical or information-theoretic properties such 
as fractal dimensions and dynamical entropies [Q]. Nevertheless, little progress has been 
made to establish a link between chaotic indicators and directly observable properties. In 
this Letter we show the existence of a general and remarkable connection: the velocity of 
correlation propagation is determined by the vanishing of the convective Lyapunov exponent. 

On the one hand, we know that chaotic systems are generally characterized by rapidly 
decaying correlations. This does not prevent the onset of sizable propagation phenomena, 
as apparent in Fig. || where, as an example, we report the space-time representation of the 
local heat-flux in a Fermi-Pasta-Ulam (FPU) chain (see later for a definition of both the 
observable and the model). A simple tool to pinpoint the existence of travelling processes 
is the autocorrelation function C{x,t) = {w(t' + t,x' + x)w{t',x')) (or, equivalently, the 
structure function) of some observable w{x, t), where x and t denote space and time variables, 
respectively, and (■) denotes a time average (we shall always assume that ergodicity holds). 
The most effective propagation processes can then be identified by estimating the velocity 
V = v that minimizes the decay of C{vt, t). In many chaotic systems the minimum is attained 
for a non-zero velocity. 

On the other hand, Lyapunov exponents represent the right tool to investigate the evo- 
lution of an infinitesimal perturbation 5{x,t). If the initial perturbation 5{x, 0) is localized 
around the origin, it has been shown that, in the limit t ^ oo, 

(5(x, t) ~ exp[A(t; = x/t)t], (1) 

where A(f) is the so-called convective or velocity-dependent Lyapunov exponent. In fact, 
it expresses the growth rate of a localized perturbation when observed in a frame moving 
with velocity v. The maximum value of A(f ) coincides with the usual maximum Lyapunov 
exponent: in convectively unstable systems, this occurs for a suitable non-zero velocity, 
while it is located in the orgin [v = 0) in spatially symmetric chaotic systems. An example 
of the typical behaviour of A(v) in this latter context can be seen in the inset of Fig. |1|. 
There, one can notice a further general feature: moving away from the maximum, A(f) 
decreases, becoming negative above some critical velocity v*. Therefore, any perturbation 
"travelling" faster than the critical velocity is exponentially damped and thus becomes 
quickly negligible. This fast damping prevents any coherent propagation of information 
and, therefore, we expect that a meaningful long-term coherence cannot be maintained for 
such large velocities. 

In the opposite limit of small velocities, it is well known that an exponential separation 
of trajectories leads to a fast decoherence, since nearby orbits rapidly enter different regions 
of the phase-space. Accordingly, strong correlations cannot again be maintained. The 
only moving frame in which neither of the two effects actively contributes to destroying 
correlations is precisely that one corresponding to the neutral point v*, so that the most 
effective propagation phenomena should occur exactly at the velocity v = v*. 

The first system where we have tested this conjecture is a chain of FPU oscillators [Q, 
an idealized microscopic Hamiltonian model for an insulating solid, 

(ji = F(gi_i - q,) - F{qi - g^+i) (2) 
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where qi represents the displacement of the ith particle from its equilibrium position and 
F{x) = —X — is the force field. This is a simple nonlinear model, introduced to test 
the ergodic hypothesis in the first numerical experiment ever performed Within the 
large number of features that have been found while investigating the dynamics of the 
above model, here we are interested in the propagation phenomena recently discovered in 
connection with the study of heat conductivity. In order to clarify the observed anomalous 
transport properties, the authors of Ref. have performed microcanonical simulations (with 
periodic boundary conditions), monitoring the local heat flux 

= - q^) + Fiq, - (3) 

From the behaviour of the autocorrelation function {j{t' + t,i' + i)j{t',i')), they clearly 
found a propagation velocity v ~ 2.47, when the energy per particle is e = 8.8. The space- 
time representation of j{t,i) in Fig. |l] demonstrates directly the existence of symmetric 
propagation phenomena along the directions identified by the cross-correlation analysis and 
denoted by the tilted arrows in the figure. 

For what concerns the spectrum of convective Lyapunov exponents, rather than letting an 
initially localized perturbation evolve, we have preferred to follow the procedure devised in 
Ref. as it suffers much less problems of finite-size corrections 0. The readers interested 
in a thorough explanation of the procedure can consult Refs. here, for the sake of 

completeness, we summarize the key steps. Very briefly, the method consists in computing 
the maximum Lyapunov exponent A(/i) of a special class of perturbations: those exhibiting 
an exponential profile with an imposed growth rate /i. This can be done by linearizing 
Eqs. (0) to obtain the standard equations for a generic perturbation 6qi. By then introducing 
the Ansatz 6qi = 6qiex]i{fii), one obtains the differential equation for the "envelope" 6qi. 
The corresponding growth rate is the generalized Lyapunov exponent A(/i). The velocity- 
dependent Lyapunov spectrum A(f ) is finally obtained by Legendre transforming A(yu) {v = 
dX/dfi; A = A(yu) — fiv). Mutatis mutandis, X{fi)/ fi can be intepreted as the "phase velocity" 
of waves with wavelength "/i" , while v can be read as the corresponding group velocity. 

The result for the FPU system is reported in the inset of Fig. |l]: it reveals two symmetric 
zeros, whose value corresponds to the direction of the straight lines visible in the corre- 
sponding pattern. However, more importantly, the absolute value of the marginal velocity 

2.46) is in full agreement with the velocity previously estimated from the behaviour of 
the correlation function. As a first check of the generality of this identity, we have repeated 
the same analysis for different energy densities, namely e = 1. and e = 0.1. The data are 
altogether reported in Table 1, where one can see that v always agrees with v* within the 
numerical error. 

The second model we have considered is the complex Ginzburg-Landau (CGL) equa- 
tion, a partial differential equation describing chaotic properties of generic systems close to 
oscillatory instabilities (it has been derived, e.g., in hydrodynamics ^ and laser physics |P), 

-9^ = ^-i^-^b)\u\\ + {l + zc)^ (4) 

where m is a complex field which, in general, represents the slowly- varying amplitude of some 
relevant mode. For b = 2.3 and c = 0, propagation phenomena are clearly visible, as it can 
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be noticed in the pattern reported in Fig. ^ (the equations have been integrated by using 



the algorithm described in Ref. |T0[ and adopting again periodic boundary conditions). It 
is instructive to notice that, at variance with the previous Hamiltonian system, now the 
pattern is no longer invariant under time reversal. This is an obvious consequence of the 
dissipative nature of the CGL equation. From the computation of the correlation function 
(|Mp(t + t', x + x')\u\'^{t', x')), we can determine the optimal propagation velocities which are 
again symmetric (see the two arrows in Fig. ^ and practically coincide with the velocity of 
the dark structures visible in the pattern. 

The comoving exponents can be again computed by linearizing Eq. (^ and assuming an 
exponential profile for the perturbation. The resulting spectrum is reported in the inset of 
Fig. ^. It is symmetric (because of the left-right symmetry in the model) and the zeros of the 
spectrum coincide with the propagation velocities of correlations (see Table 1) represented 
by the two tilted arrows in the same figure. 

The last system where we have tested our conjecture, is an experimental one, namely, a 
CO2 laser with delayed feedback. This is a physical system that has been used in the past to 
investigate several instabilities both in the regime of low and high-dimensional chaos |Tl],|12[. 
Although this is not, strictly speaking, a spatially extended system, it can be interpreted 
as such. The idea consists in decomposing the time variable a.s t = n + st |]TB|], where the 
parameter r is close to the delay time (fixed equal to 400 fis in the experiment), n = Int(t/r) 
is the new, discrete, time variable, and s (0 < s < r) is a space-like variable. In a sense, 
this system is complementary to a chain of oscillators, since the discrete and continuous 
character of time and space axes are exchanged. The validity of this decomposition has 



been established in Ref. [0]; here, the reader can appreciate its meaningfulness by looking 
at the space-time representation of the intensity of the emitted radiation in Fig. |^. While 
looking at the pattern, it is important to realize that the slope of the various coherent 
structures depends on the value of r adopted in the time decomposition. Here, for the 
sake of clarity, we have adjusted r in such a way that the pseudo-rolls, responsible for the 
optimal correlations, are almost vertical (see arrow 1). A further important difference with 
the previous models is the absence of left-right symmetry. Even more, causality implies that 
no leftwards (negative) propagation is possible at all. 

Since no theoretical model exists which reproduces the laser djTiamics with a sufficient 
accuracy, one must exclusively rely on the experimental data. However, it has been recently 
developed an approach to reconstruct the dynamics of delayed systems, even when it is so 
high-dimensional that the standard embedding technique is bounded to fail |TB|,|TB[. In short. 



given a time-series ?/„, the method consists in constructing an embedding space composed 
of two-window vectors, 

Vn = iVn, Vn-l, ■ ■ ■ , Vn-m+l, Un-Ti ■ ■ ■ i Vn-T-m+l), (5) 

where the distance T between the two windows coincides with the delay, while the length 
m is the effective number of variables involved in the dynamics. Once the proper values 
of both T and m have been identified (by minimizing the forecast error), a local linear 
model can be constructed, by fitting the behaviour in a suitable neighbourhood of each 
point in the embedding space. In the case of the CO2 laser, m = 5 guarantees an excellent 
reproduction of the original dynamics ||T^. From the knowledge of the empirical model, 
one can again compute the convective Lyapunov exponents by going through the same 
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intermediate steps. The corresponding spectrum is reported in the inset of Fig. | (notice 
that, given the pecuharity of the space-time reconstruction, velocities are here expressed in 
time units - more precisely us per number of delay units), where we see that it is restricted 
to the positive-w domain (see Ref. fl^). As a consequence of this asymmetry, both zeros 
are positive. Thus, in principle, one might expect two different propagation velocities. The 
correlation analysis of the experimental data has instead revealed a single velocity, which 
coincides with the smallest zero of the convective Lyapunov spectrum (see Table 1). We 
cannot rule out the possibility that a second, much less effective, propagation exists with 
a larger velocity in coincidence with the second zero. However, even if this is not the 
case, we can at least state that whenever propagation of correlations can be inferred from 
the evolution of some observable, it must coincide with a netrually stable velocity for the 
spectrum of convective Lyapunov exponents. What are the (possibly) additional ingredients 
ensuring the actual existence of propagation phenomena and determining their strength 
seems to be a much harder problem that we plan to attack in the future. 

We conclude, by restating that we have found a clear evidence of a strict link between 
a "large" -scale feature like the optimal propagation of correlations and the evolution of 
infinitesimal perturbations, a problem which, because of its very nature, can be formulated 
in terms of Lyapunov exponents. Accordingly, albeit deterministic chaos contributes to a 
fast decay of correlations in low-dimensional systems, it is compatible with propagation 
phenomena as soon as spatial degrees of freedom come into play. 
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TABLES 



System 


V* 


V 


FPUa 


2.46 ±0.01 


2.47 ±0.04 


FPUb 


1.55 ±0.01 


1.53 ± 0.05 


FPUc 


1.11 ±0.01 


1.2 ±0.15 


CGL 


0.38 ±0.01 


0.37 ±0.02 


Laser 


(4.9 ± 0.1)/x sec 


(4.8 ± 0.2)/Lt sec 



TABLE L The velocity v* of the zero convective Lyapunov exponent and the velocity v of the 
optimal propagation of correlations in an FPU chain (for e = 8.8 (a), e = 1 (b), and e = 0.1 (c)), 
in the CGL equation, and in the CO2 laser system (from experimental data). Except for the last 
case, velocities are expressed in adimensional units. 
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FIGURES 




FIG. 1. Space-time representation of the heat-flux dynamics in a chain of 512 FPU oscilla- 
tors. Time flows from bottom to top as in the other figures. The total time span is 51.2 in the 
rescaled units of Eq. (1). The two tilted arrows correspond to the optimal propagation velocity 
as determined from the computation of the space-time correlation function of the heat-flux (see 
Eq. H). The vertical arrow corresponds to the direction of the maximal growth rate for infinitesimal 
perturbations. In the inset, the dependence of the convective Lyapunov exponent is reported as a 
function of the velocity, i.e. the slope in the space-time representation. 
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FIG. 2. Space-time representation of the square amplitude of the complex Ginzburg-Landau 
equation in a system of length 51.2 over a time span equal to 50. The meaning of the arrows and 
of the inset is the same as in the previous figure. 



9 




FIG. 3. Space-time representation of the field intensity in the CO2 laser experiment described 
in the text. The horizontal axis represents the time variable within a delay unit of length 406 /isec, 
while the y direction corresponds to different units (see the text for a more detailed explanation). 
The central arrow again denotes the direction characterized by the maximal growth rate. 
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